Reducing the Linewidth of an Atom Laser by Feedback 
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A continuous atom laser will almost certainly have a linewidth dominated by the effect of the 
atomic interaction energy, which turns fluctuations in the condensate atom number into fluctuations 
in the condensate frequency. These correlated fluctuations mean that information about the atom 
number could be used to reduce the frequency fluctuations, by controlling a spatially uniform poten- 
tial. We show that feedback based on a physically reasonable quantum non-demolition measurement 
of the atom number of the condensate in situ can reduce the linewidth enormously. 
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An atom laser is a continuous source of coherent atom 
waves, analogous to an ordinary laser which is a con- 
tinuous source of coherent photon waves (light) jl],|]]. 
Ideas for creating an atom laser were published indepen- 
dently by a number of authors &-@], shortly after the 
first achievement of Bose-Einstein condensation (BEC) 
of gaseous atoms |t]-|9| . There have since been some im- 
portant advances in the coherent release of pulses [ ^0|JTl| 
and beams Jl2|,[l3| of atoms from a condensate. Since the 
condensate is not replenished in these experiments, the 
output coupling cannot continue indefinitely [p3. Nev- 
ertheless they represent major steps towards achieving a 
continuously operating atom laser. 

The coherence of an atom laser beam can be defined 
analogously to that of an optical laser beam: the atoms 
should have a relatively small longitudinal momentum 
spread, they should ideally be restricted to a single trans- 
verse mode and internal state, and their flux should be 
relatively constant (|] . A fourth condition, rarely consid- 
ered for optical lasers because it is so easily satisfied, is 
that the laser beam be Bose degenerate. This requires 
that the atom flux be much larger than the linewidth 
(the reciprocal of the coherence time) . A crucial con- 
tributor to the linewidth of an atom laser is the colli- 
sional interaction of atoms, which is negligible for pho- 
tons. This is difficult to avoid because it is the collisions 
between atoms that enables BEC by evaporative cooling, 
at present the only method for achieving an atom laser. 

In this Letter we show that using a feedback mecha- 
nism can reduce the effect of atomic interactions on the 
atom laser linewidth by a factor as large as the square 
root of the atom number. For a single-mode condensate 
the dominant effect of collisions is to turn atom number 
fluctuations in the condensate into fluctuations in the en- 
ergy, which are equivalent to frequency fluctuations. By 
monitoring the number fluctuations, it is possible using 
feedback to largely compensate for the linewidth caused 
by these frequency fluctuations. The key practical points 
are that the measurement does not rely upon any exter- 
nal condensate phase reference, and that the control re- 
quires only the ability to change the energy of the atoms, 



which could be done with a spatially uniform optical or 
magnetic field. 

We begin by deriving the standard laser linewidth (for 
non-interacting bosons) using a simple method which is 
applied to all later cases. We then derive the broad- 
ened linewidth for an atom laser with strongly interact- 
ing atoms. Finally, we show that feedback based on a 
quantum non-demolition (QND) measurement of atom 
number in the the condensate can greatly mitigate this 
linewidth broadening. 

(a) Standard Laser Linewidth. To derive the standard 
linewidth we use the usual single- mode model of the laser 
p5|^6| . Far above threshold, the laser mode has Poisso- 
nian number statistics. In the absence of thermal or other 
excess noise, its dynamics are modeled by the completely 
positive master equation ]l7|,|il| 



p = K^P[a t ]^l[a t ]"V + KV[a]p. 



(1) 



Here n is the loss rate, /i» 1 the stationary mean bo- 
son number, and a the annihilation operator for the laser 
mode. The superoperators V and A are defined as usual: 



V[A]B = ABA^ - A[A}B, 
A[A]B = (A^AB + BA^A)/2. 



(2) 
(3) 



That Eq. (fil) is of the Lindblad form follows from the 
identity |,0 £>[a t ].4[ at ]~ 1 = J °° dqD[ate- qaai / 2 ]. The 
stationary solution to this master equation is 
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where r — ^fp and \re l8 ) is a coherent state jl5| , |i~f| . 

The coherence time of a laser is roughly the time for 
the phase of the field to become uncorrelated with its ini- 
tial value. It is determined by the stationary first-order 
coherence function 



.gW(i) = Tr[a t e £t ap ss ]/Tr[p ss a t a] 
A simple and useful definition is B 
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Here the 1/2 is so that, for the standard laser, I = r ( 
will be the standard linewidth. In the cases we consider 
it is a very good approximation jl9] to put Ifif^^t)! = 
e -iwotg(i)^^ f or sorne frequency ljq. From Eq. (||), this 
allows the integral in Eq. (0) to be evaluated, yielding 



Tcoh — -Tr[a t (£ - iuj ) 1 ap ss ]/2Tr[p ss a t a]. 



(7) 



Equation (Q) can be evaluated numerically for ex- 
ample using the Matlab quantum optics toolbox poj . 
Analytically it is easier to use the fact that Eq. (p3) 
is unchanged if p ss is replaced by the coherent state 
\re l9 ){re l9 \ for arbitrary 9. Using any suitable phase- 
space (a, a*) representation, the coherence function then 
becomes g^(t) = (a*(t)) / (a*(0)), where a(0) has a dis- 
tribution corresponding to |re* e ). Because fluctuations in 
the intensity n = \a\ 2 are relatively small in a laser with 
/i>l, the coherence time is well approximated by 
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Here 8V$(t) is the variance of Scf)(t) — arg[a(t)/a(0)], and 
the second approximation relies on 5<fi(t) having Gaussian 
statistics, as will be justified below. 

For our laser model, the Q-function is the most conve- 
nient representation because of the identity 



tl-i, 



t(- 



k=l 



Q{n,4>), 



(9) 



which, since the higher order derivatives are negligible, 
can be truncated at k — 2. The master equation ([!]) thus 
turns into a Fokker-Planck equation (FPE) for Q(n, (f>) 
which can be linearized. Under this FPE, the number 
statistics remain those of the initial coherent state, as 
does the number-phase covariance C n <p (he. it remains 
zero), and the phase has Gaussian statistics with a vari- 
ances that increases as 8V<j,(t) — («/2^)t. Substituting 
this into Eq. (@), we obtain r = 2/i/k. This time is pre- 
cisely the coherence time. Its reciprocal is the standard 
laser linewidth £ = rc/2/z |J|-[l||. 

(b) Atom Laser Linewidth. As a model for an atom 
laser we take the standard laser master equation (|l|) and 
add a term —iC{a) a' aa, p\. This represents the self- 
energy of the atoms due to collisions, where 



( , = 2nna 1 j d 3 r |^ (r) |4 



(10) 



Here ip{ r ) is the wavefunction for the condensate mode, 
and a s is the s-wave scattering length. For the experi- 
ments in Refs. |]To|— |l3| , ip(r) can be determined using the 
Thomas-Fermi approximation, and we use this to obtain 



the numerical values which appear later. The Hamilto- 
nian CoJo) aa has no effect on the atom number statistics. 
Linearizing the resultant FPE for the Q function yields 
the phase-related second-order moments 



C n ^) = |(e- rei -l), 
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(11) 
(12) 



Here, we have used x = 4/xC / n as a dimensionless pa- 
rameter for the strength of the atomic interactions. 

This expression for 8V$(t) implies that the integrand 
g^{t) in Eq. (||) has the same structure as the anal- 
ogous expression, Eq. (184), derived in Ref. [|l). This 
was for a condensate in dynamical equilibrium with ther- 
mal atoms. Moreover, the three time scales identified 
m Ref. H have the same physical origins as those in 
Eq. (|l^). The integral in Eq. (g) can be evaluated ana- 
lytically in two limits: 
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K(l + x 2 )/2fi for X< 
2kx/\ / 2ttJ1 for x > \fP- 
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These correspond to the characteristic time constants 
of Eqs. (187) and (186) of Ref. @ respectively (see 
Ref. p^] for a further discussion). Our two expressions 
agree at x — \/8^/7r, so we have an expression for how 
I scales for all values of x- The second expression for 
t, in the regime where the nonlincarity is dominant, is 
familiar as the collapse time of an initial coherent state 
in the absence of pumping or damping ]23| , p4| ]. 

Using the preliminary atom laser experiments |l0|-|l3|] 
as a guide to realistic parameter values, the dimensionless 
interaction strength x l& found to always satisfy x 3> 1. 
This implies a linewidth for the atom laser far above the 
standard limit. If x ^ A^ 2 , then I would be larger than 
the output flux K/i, and the output would cease to be 
coherent. It is thus of great interest to find methods for 
reducing the linewidth due to atomic interactions. 

(c) Effect of QND-based Feedback. Atomic interac- 
tions do not directly cause phase diffusion. Rather, they 
cause a shearing of the field in phase space, with higher 
amplitude fields having higher energy and hence rotat- 
ing faster. The linewidth-broadening which results is a 
known effect for optical lasers with a Kerr (x^) medium 
p5| . The shearing of the field is manifest in the finite 
value acquired by the covariance C n< j,(t) in Eq. (|TT|) . It 
means that information about the condensate number is 
also information about the condensate phase. Hence, we 
can expect that feedback based on atom number mea- 
surements could enable the phase dynamics to be con- 
trolled, and the linewidth reduced. 

It might be thought that one could measure the atom 
number by diverting some of the atom laser output beam. 
It turns out that this sort of measurement is effectively 
useless for reducing the linewidth. For this reason we 
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consider instead quantum non-demolition (QND) mea- 
surements of atom number, which works very well |26] ]. 

QND atom number measurements can be performed 
on the condensate in situ using dispersive imaging tech- 
niques |2|] . We consider a far-detuned probe laser beam 
of frequency uj p and cross-sectional area A (assumed 
larger than the condensate) which passes through the 
condensate. For simplicity, we will assume that the dis- 
tortion of the beam front, and the mean phase shift, are 
removed by a suitable "anti-mean-condensate" lens. For 
a single atom, the interaction Hamiltonian with the probe 
beam in the large detuning limit is 
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(14) 



where A, T and I sa ± have their usual meaning. Here 
b is the annihilation operator for the probe beam, nor- 
malized so that Tkjpb'b is the beam power. The effective 
interaction Hamiltonian for the whole condensate, minus 
the mean phase shift, is thus 



K6(ala-(ipb, 



(15) 



where 8, defined in Eq. (14), is the phase shift of the 
probe field due to a single atom. 

From Eq. (|l5[), the back-action of the probe fluctua- 
tions on the condensate can be evaluated using the tech- 
niques of Ref. |E9[ , and results in the extra phase diffusion 



(3 2 V[e- l9 ^ a ~^]p ~ NV[a)a]p. 



(16) 



Here the approximation requires ^fpd -C 1, and 

N = 2 9 2 = P0 2 /hLu p , (17) 

where (3 = {b- ln ) and P is the mean beam power. Equa- 
tion (|lj) also gives the output probe field |2{| 



= e- ie ^ a -"h in ~ b- m - 10(0*0 - ft)p, 



(18) 



where the same approximation has been used. The con- 
densate number fluctuations can thus be measured by 
homodyne detection of the output phase quadrature 



= -i&out + ibi 



2(39{a)a-p). (19) 



In order to control the phase dynamics of the conden- 
sate, we wish to use the homodyne current to modulate 
the phase. This can be done, for example, by applying 
a uniform magnetic field or far-detuned light field across 
the whole condensate. We model this by the Hamiltonian 



H ib (t) = -ha)aJ — h(s)lj ui (t-s)ds, (20) 
V Jo 

where h(s) is the response function of the feedback 
loop and A is a dimensionless measure of the feedback 



strength. Neither the measurement nor the feedback af- 
fect the atom number statistics. In the ideal limit of in- 
stantaneous feedback, h(s) — > S(s) and we can apply the 
Markovian theory of Ref. . The total master equation 
for the atom laser is then 



p = K//D[o t ].A[a t ]~V 



kD[cl\p — iC[i 



,t„t, 



+ iiVy^aVaa, p] + N ^1 + -^j V[a ] a]p. (21) 

Here we have allowed for a detection efficiency rj [p9| , 
and dropped terms corresponding to a frequency shift, 
and defined a dimensionless parameter v = ApN/n. 
Proceeding as before, we find the Q phase variance: 
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(22) 



If the feedback is to reduce the linewidth, we need 
(x — \f~Xv) 2 <C 4/i and rj must be large enough that 
tj/i » 1. In this case, Eq. (||) has a simple analytical 
solution: 



1 = 



K 

4^ 



(x-v^a)' 



(23) 



Assume x is large enough that 2^/rj\ > 1 , which is the 
physically interesting regime where the self-energy is im- 
portant. Then £ is minimized for a feedback strength of 
A = y^x— 1/2, and a measurement strength of v = A/77. 
The optimum linewidth in this case is 



2ft 



VV 4/7 



(24) 



Note that for a large interaction strength \ there is an op- 
timum measurement strength v ~ x/\/v (or equivalently 
N ~ C/^/rj) independent of the output coupling rate k. 
From Eq. (Ell) we see that this effectively cancels the self- 
energy term in the master equation (since A = r\v) . A 
weak measurement will give poor information about the 
atom number, with a high noise-to-signal ratio. Feeding 
back the noisy current to counter the shearing caused by 
the nonlinearity will thus add large phase fluctuations to 
the condensate. On the other hand, if the measurement 
is too strong the measurement back action in the form 
of phase diffusion, as discussed above Eq. ([Tt]), will itself 
dominate the linewidth. 

In Fig. 1 we plot the approximate analytical expres- 
sions for the linewidth in the absence [Eq. (|l3|)1 and pres- 
ence [Eq. (pi])] of feedback as a function of nonlinearity 
X- We also plot exact numerical results using Eq. ©, for 
p = 60. The agreement is very good. It is evident that 
the QND feedback offers a linewidth much smaller than 
that without feedback for all values of x > l/2y/V- In fact 
for x> 1 the reduction factor approaches a maximum of 
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y/8fj./iT. Most importantly, with feedback, the laser out- 
put remains coherent until X ~ M 2 ; a much higher value 
than the /x 3 / 2 which applies in the absence of feedback. 




X 

FIG. 1. Atom laser linewidth for k = n = 1 and /i = 60, 
plotted with and without feedback using both analytical and 
numerical methods (see text). 

Let us summarize. An atom laser will almost certainly 
suffer great linewidth broadening due to the collisional 
self-energy of the atoms. This is because the self-energy 
produces a correlation between the number and phase 
fluctuations of the condensate. We have shown that this 
broadening can be enormously reduced by controlling the 
phase of the condensate based on a QND measurement 
of the number of atoms in the condensate. The factor 
of reduction can be as large as the square root of the 
number of atoms in the condensate (that is, a factor of 
perhaps 10 3 ). 

A question of interest is, how easy is it to obtain a QND 
measurement of sufficient strength to optimize the feed- 
back? We have seen above that we require v ~ x, which 
is equivalent to N ~ C. For the typical BEC parameters 
of Refs. C - lfrV 1 . It may be verified from 

Eq. (|l7| ) that it is very easy to obtain a measurement 
strength N this large, even with A 3> T. 

A related question is, how much of a problem is atom 
loss due to spontaneous emission by atoms excited by the 
detuned probe beam? The rate of this loss (ignoring re- 
absorption) is r x (Excited Population). We would like 
the ratio of this loss rate to the output loss rate K[i to be 
small. In the A ^> T limit, this ratio is given by 

r (r/2A) 2 (//2/ sat y ^ 2 x I sat A 

For physically reasonable parameters of x ^ 10 3 , 7 sat ~ 
10W/m 2 , A ~ I0- n m 2 , ol ~ 3 x 10 15 s - \ T - 3 x 
10 6 s-\ and \i ~ 10 6 , Eq. @ is indeed small (~ HT 1 ). 

In conclusion, there appear to be no fundamental rea- 
sons that this proposal could not be put to good use when 



a continuous atom laser is realized. 
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